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Abstract 



■^ i In this paper, the bending and the free flexural vibration behaviour 

7 j I of sandwich functionally graded material (FGM) plates are investi- 

(-H ■ gated using QUAD-8 shear flexible element developed based on higher 

order structural theory. This theory accounts for the realistic variation 
of the displacements through the thickness. The governing equations 
obtained here are solved for static analysis considering two types of 
sandwich FGM plates, viz., homogeneous face sheets with FGM core 
h>. I and FGM face sheets with homogeneous hard core. The in-plane and 

■^ ' rotary inertia terms are considered for vibration studies. The accu- 

[y^ . racy of the present formulation is tested considering the problems for 

^SJ ! which three-dimensional elasticity solutions are available. A detailed 

^f^ ' numerical study is carried out based on various higher-order models 

^^ ■ to examine the influence of the gradient index and the plate aspect 

^^ ! ratio on the global/local response of different sandwich FGM plates. 

Keywords: Functionally graded material plate, higher-order theory, thermal loading, 
^^ , mechanical loading, free vibration, shear flexible element. 

1 Introduction 

With ever- increasing demand for high strength-to-weight ratio materials, the 
engineered materials have been replacing the conventional materials in auto- 
motive, nuclear and aerospace industries. In general, most of the engineered 
materials are inspired from nature. A new class of materials was introduced 
by Japan scientists [1] to decrease the thermal stresses in the propulsion 
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systems and the airframe of the space planes. This class of engineered ma- 
terials was coined as functionally graded material (FGM). These materials 
are made up of mixture of ceramics and metals, that are characterized by 
the smooth and continuous variation in the properties from one surface to 
another [H |2] . For the structural integrity, the FGMs are preferred over the 
fiber-matrix composites that may results in debonding due to the mismatch 
in the mechanical properties across the interface of two discrete materials 
bonded together. This has attracted lot of researchers in understanding the 
mechanics and mechanism of the FGM structures. 

Among the various structural constructions, the sandwich types of struc- 
tures are commonly used in the aerospace vehicles, because of its outstand- 
ing bending rigidity, low specific weight, excellent vibration characteristics 
and good fatigue properties. Laminated composite types of constructions 
are, in general, adopted in sandwich structures. However, due to the sudden 
change in the material properties from one layer to another, the variation of 
the interfacial stress distribution is significant at the facesheet-core interface. 
Furthermore, the response of such laminated composites depend on the bond- 
ing characteristics. In contrast, the FGM sandwich can alleviate the large 
variation in the interfacial stress distribution, because of the gradual varia- 
tion of the material properties at the facesheet-core interface. For predicting 
the realistic structural behaviour of such components, more accurate analyt- 
ical/numerical analysis based on the three-dimensional models may be com- 
putationally involved and expensive. Hence, among the researchers, there is 
a growing appreciation of the importance of applying two-dimensional theo- 
ries with new kinematics for the evolution of the accurate structural analysis. 
Various structural theories proposed for the FGM structures have been ex- 
amined and some of the important contributions pertaining to the sandwich 
FGM plates are discussed here. The most commonly used FGM sandwich 
constructions are: the FGM facesheet with homogeneous core and the ho- 
mogeneous facesheet with FGM core. These sandwich constructions can be 
considered for the requirement of light weight and high bending stiffness in 
design (for example, the placement of actuators/sensors in the field control 
and thermal/mechanical load bearing component design) by appropriately 
selecting the soft/hard core metal or the ceramic layers. 
Venkataraman and Sankar p] and Anderson |4j have studied the effect of 
FGM core in sandwich beam on the interfacial shear stresses. Pan and 
Han [5] analysed the static response of the multilayered rectangular plate 
made of functionally graded, anisotropic and linear magneto-electro-elastic 



materials. Das et al., [B] have investigated a sandwich consisting of FGM soft 
core with relatively stiff orthotropic facesheets employing triangular plate el- 
ement. Ganesan et al, \J] studied static and dynamic response of FGM plate 
with viscoelastic core, whereas Shen |8] examined FGM sandwich plates with 
piezo-electric core subjected to thermo-electro-mechanical loading. Zenk- 
our [HI [ID] studied analytically the static and dynamics of sandwich FGM 
plates with homogeneous hard core using sinusoidal shear deformation the- 
ory. Li et al. [11] have presented three-dimensional analytical solutions for 
multi-layer FGM sandwich plates based on Ritz method in conjunction with 
Chebhyshev polynomial series. It is observed from these studies that first- 
order, third-order, and sinusoidal shear deformation theories have been ex- 
tensively used for the analysis of sandwich FGM plates. However to author's 
knowledge, the theories accounting the variation of in-plane displacement 
through the thickness, and the possible discontinuity in slope at the inter- 
face, and the thickness stretch affecting transverse deflection is not exploited 
while investigating the structural behaviour of FGM sandwich structures. A 
Layer wise theory is the possible candidature for this purpose, but it may be 
computationally expensive as the number of unknowns to be solved increases 
with the increase in the number of mathematical or physical layers. Ali et 
al. [12], and Ganapathi and Makhecha [13] have employed a new higher-order 
plate theory based on global approach, for multi-layered laminated compos- 
ites by incorporating the realistic through the thickness approximations of 
the in-plane and transverse displacements by adding a zig-zag functions and 
higher-order terms, respectively. This formulation has proved to give very 
accurate results for the composite laminates. Such model for the current 
problems is now explored as a candidature while comparing with the three- 
dimensional model. 

In this paper, a C° 8-noded quadrilateral plate element with 13 degrees of 
freedom per node [T3| [H] based on the higher order theory [12j is employed 
to study the static deflection and the free vibration analysis of thick/thin 
sandwich functionally graded material plates. The efficacy of the present for- 
mulation, for the static analyses subjected to the thermal/mechanical loads 
and the free vibration analysis is illustrated through the numerical studies. 
The accuracy of the present element with those of the other two-dimensional 
numerical/analytical models and the three-dimensional elasticity analysis are 
discussed considering various parameters. 

The paper is organized as follows. The next section will give an introduction 
to the FGM and a brief overview of the higher order theory. Section §3] 



describes the 8-noded quadrilateral plate element used in the current study. 
Section §1] present numerical results for the static deflection and the free 
vibration of thick/thin sandwich FGM plates, followed by concluding remarks 
in Section 851 



2 Theoretical Formulation 

2.1 Functionally graded material plate 

Consider a rectangular FGM plate with co-ordinates x, y along the in-plane 
directions and z along the thickness direction as shown in the Figure [H The 
material is assumed to be graded only in the thickness direction according to a 
power-law distribution. The homogenized material properties are computed 
using the rule of mixtures. 




Figure 1: Co-ordinate system of rectangular FGM plate. x,y along the in- 
plane directions and z along the thickness direction. 



Estimation of mechanical properties 

The effective material properties for each layer, viz.. Young's modulus, Pois- 
son's ratio and the mass density is estimated by the following power-law 
function: 



P{Z) = PrnVm + PcVc (1) 

Here, Vi {i = c, m) is the volume fraction of the phase materiaL The sub- 
scripts c and m refer to the ceramic and the metal phases, respectively. Note 
that Pm and Pc are the properties of the metallic and the ceramic phases, 
respectively. The volume fractions of the ceramic and the metal phases are 
related by Vc + Vm = I, and Vc is expressed as: 

where n in Equation (|2]) is the volume fraction exponent, also referred to as 
the gradient index in the literature. The properties of the FGM plate vary 
continuously through the thickness based on a power-law function, given 
by Equation ([T]). In this study, the following two types of power-law FGMs 
are considered. 

Type A - FGM face sheet and homogeneous hard core 

In this case, the volume fraction of the FGMs is assumed to follow 
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Vs{z)=[ ^ , ze[z,,z4 (3) 

\Z3 — Z4/ 

where Vi{i = 1,2,3) is the volume fraction of layer i and n is the gradient 
index. The core is considered as a fully ceramic material. The top and 
bottom surfaces of the plate are metal-rich. 

Type B - Homogeneous face sheet and FGM core 

In this case, the volume fraction of the FGMs is assumed to follow 
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The core is considered as a functionally graded material. The top surface of 
the plate is metal rich and the bottom surface of the plate is ceramic rich. 
Figure H] shows the variation of the volume fractions of ceramic and metal, 
respectively, in the thickness direction z for the two types of sandwich FGM 
plates considered in this study. 

2.2 Plate formulation 

The sandwich FGM plate is assumed to be made up of three discrete layers. 
The material properties for fc*'' layer is governed by Equation ([3]) or Equation 
(jlj). The in-plane displacements u^ and v'', and the transverse displacement 
w^ for the k^^ layer, are assumed as 



u''{x,y,z,t) = Uo{x,y,z,t) + ze^^x.y.t) + z'^(3r,{x,y,t) + 

z^(t)x{x, y, t) + S^ip^cix, y, t) 

v''{x, y, z, t) = Voix, y, z, t) + z9y{x, y, t) + z^Py{x, y, t) + 

z^(j)y{x,y,t) + S^iljyi^x.y.t) 
w^{x, y, z, t) = Woix, y, t) + zwi{x, y, t) + z'^T{x, y, t) (5) 

The terms with even powers of z in the in-plane displacements and odd pow- 
ers of z occurring in the expansion for w^ correspond to the stretching prob- 
lem. However, the terms with odd powers of z in the in-plane displacements 
and the even ones in the expression for w'' represent the flexure problem. 
Uo, Vo and Wo are the displacements of a generic point on the reference sur- 
face; 6x and 6y are the rotations of the normal to the reference surface about 
the y and x axes, respectively; Wi, (3^, f3y,T,(j)x and (j)y are the higher-order 
terms in the Taylor's series expansions, defined at the reference surface. Vx 
and ipy are generalized variables associated with the zigzag function, S^. The 
zigzag function, S^, as given in [151 [ISl US], is defined by 
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Figure 2: Through thickness variation of volume fraction for various types 
of sandwich FGM plates for gradient index n =2. 



gk ^ 2(-l)^^ (6) 

where Zk is the local transverse coordinate with the origin at the center of 
the A;*'* layer and h/^ is the corresponding layer thickness. Thus, the zigzag 
function is piecewise linear with values of -1 and 1 alternatively at different 
interfaces. The 'zig-zag' function, as defined above, takes care of the inclusion 
of the slope discontinuities of u and v at the interfaces of the sandwich 
plate as observed in the exact three-dimensional elasticity solutions of thick 
sandwich functionally graded materials. The main advantage of using such 
formulation is that such a function is more economical than a discrete layer 
approach [T71 [18] . 

The strains in terms of mid-plane deformation, rotations of normal and higher 
order terms associated with displacements are: 

t } (7) 

The vector £bm includes the bending and membrane terms of the strain com- 
ponents and vector Ss contains the transverse shear strain terms. These 
strain vectors are defined as [121 [ISI [Ti 
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and. 



The subscript comma denotes partial derivatives with respect to the spatial 
coordinate succeeding it. The constitutive relations for an arbitrary layer k 
can be expressed as: 



"" \ ^xx ^yy ^ zz ^xy ^xz ^yz J 

= Q*^! £bm Ss f (12) 

where Qk is the stiffness coefficient matrix defined as: 



E(z)_ , _ v{z)E{z) 



Vll — V22 ~ 1 9' V12 ~ T / \9' ^16 ~ V26 ~ '-' 



Q44 - Q55 - ^66 - 2(1 + ^(^)) (13) 



where the modulus of elasticity E{z) and Poisson's ratio ulz) are given by 
Equation ([T]). The governing equations of motion are obtained by applying 
Lagrangian equations of motion given by: 
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where 6i is the vector of degrees of freedom and T is the kinetic energy of 
the plate given by: 
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dz 
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where pk is the mass density of the k^^ layer, hk and hk+i are the z coordinates 
corresponding to the bottom and top surfaces of the /c*'* layer. The total 
potential energy function U is given by: 



U{S) = - 



22 / ^^^ dz 



dxdy 



qw dxdy 



(16) 



where q is the distributed force acting on the top surface of the plate. Sub- 
stituting Equations ( IT6l) and ( TT5l) in Equation (IHj), one obtains the following 
governing equations for static deflection and free vibration of plate. 



Static deflection 



Kd = i 



(17) 



Free vibration 

M5 + K5 = (18) 

where M is the mass matrix, K is the stiffness matrix and f is the external 
force vector. In the present study, while performing the integration, terms 
having thickness co-ordinate z are integrated with higher order Gaussian 
quadrature, because the material properties vary continuously through the 
thickness. The terms containing x and y are evaluated using full integra- 
tion with 3 x3 Gauss integration rule. The frequencies and mode shapes 
are obtained from Equation flTSl) using the standard generalized eigenvalue 
algorithm. 
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3 Element description 

In this paper, C° continuous, eight-noded serendipity quadrilateral shear flex- 
ible plate element with 13 nodal degrees of freedom 

{uo, Vo, Wo, Ox, 6y, wi, (3x, Py, r, (px, (j)y, ipx, ipy : 13 DOF is used). The finite ele- 
ment represented as per the kinematics based on Equation (I5j) is referred to 
as Q8-HSDT13 with cubic variation. Interested readers are referred to the 
literature [131 113] ? where the behaviour of the element is described in detail. 
The element is shown to be free from locking syndrome, absence of spuri- 
ous energy modes, passes patch test and exhibits faster convergence p^lT^. 
Three more alternate discrete models are proposed to study the influence of 
higher-order terms in the displacement functions, whose displacement fields 
are deduced from the original element by deleting the appropriate degrees of 
freedom. These alternate models, and the corresponding degrees of freedom 
are listed in Table [1] 

Table 1: Alternate eight-noded finite element models 
Finite element model Degrees of freedom per node 

HSDT13 (present) Uo, Vo, Wo, 6^, By, wi, (5^, Py, T, (p^,, (py, ipx, ipy 

HSDTll Uo, Vo, Wo, Ox, Oy, Wi,[5x, Py, r, (px, (Py 

HSDT9 Uo,Vo,Wo,0x,0y,l3x,l3y,(px,(Py 

FSDT Uo,Vo,Wo,9x,9y 



4 Numerical results and discussion 

In this section, we present the static response and the natural frequencies 
of sandwich FGM plates using the eight-noded shear flexible quadrilateral 
element. The effect of plate aspect ratio and gradient index are studied. In 
this study, only simply supported boundary conditions are considered and 
are as follows: 



Uo = Wo = 9x = Wi = T = (3x = (px = -ipx = 0, 


on y = 0,b 


Vo = Wo = Oy = Wl = T = f3y = (py = ^Py = 0, 


on X = 0,a 



(19) 
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where a and h refer to the length and width of the plate, respectively. The 
FGM plate considered here consists of Alumina and Aluminum. The mass 
density p and Young's modulus E are: pc = 3800 kg/m^, Ec = 380 GPa for 
Alumina and pm = 2707 kg/m^, Em = 70 GPa for Aluminum. Poisson's ratio 
z/ is assumed to be constant and taken as 0.3 for the current study. Here, 
the modified shear correction factor obtained based on energy equivalence 
principle [T^l ED] is used for the FSDT model. The transverse shear stresses 
are evaluated by integrating the three-dimensional equilibrium equations for 
all types of elements. For the current study, three different core thickness 
(1-1-1, 1-2-1, 2-2-1) for both Type A and Type B sandwich FGM plates, 
three thickness ratio a/h (5,10,100) and four gradient indices n (0,0.5,1,5) 
are considered. 

4.1 Static analysis 

The static analysis is conducted for Type A FGM sandwich plate. The 
following two types of loading are considered: 

• Mechanical loading: q = qoSin{7Tx/a)sin{7Ty/b), 

• Thermal loading: T = To{2z/h) sin(7rx/a) sini^ny/b). 

where Qo and To are the amplitudes of mechanical and thermal loads, respec- 
tively. The physical quantities are nondimensionalized by relations: 



lOO^o, , 



\(^xx: ^yyi ^xy) 



QohS^ 

_ lOOEo, . 

y^xxi '-'yyi ^xy) 



qoS' 

XZ1 

QoS 



{^x..^y.) = ^'''':y (20) 
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for the applied mechanical load and by: 



hamioo 
w = ^ (x) 



Kyxxi ^yyi ^xy) 



\^xxi ^yyi ^i 



yy, "xyj 



^m^m-'- o 



{Oxz, Oyz) = p ^ (21) 

for the applied thermal load, where the reference value is taken as Eq = 1 
GPa, S = a/h and Em, ctm are the Young's modulus and the co-efficient of 
thermal expansion corresponding to the metallic phase. Based on a progres- 
sive mesh-refinement, an 8 x 8 mesh is found to be adequate to model the 
full FGM plate for the present analysis. Before proceeding to the detailed 
analysis of the static response of the FGM sandwich plate for the applied me- 
chanical and thermal loads, the present formulation is validated considering 
the problems for which three-dimensional elasticity solutions are available. 
In Table [2] the converged displacements and stresses obtained for the Al/SiC 
functionally graded square plate under mechanical loading is compared with 
the three-dimensional elasticity solutions |2T]. The effective material prop- 
erties are based on the Mori-Tanaka homogenization scheme. The material 
properties for Al are: Em = 70 GPa, Um = 0.3, am = 23.4 xlO~^ /K and 
Km = 233 W/mK. For SiC, the material properties are: Em = 427 GPa, 
Um = 0.17, am = 4.3 xlO~^ /K and Km = 65 W/mK. It can be seen that 
the results from the present formulation are in good agreement with the 3D 
elasticity solutions. 

Next, the numerical study is carried out for Type A FGM sandwich plate 
under mechanical/thermal loads for two gradient indices, n (= 0,0.5) and for 
different thickness of the plate. Numerical results are tabulated in Tables [3] 
- m for an applied mechanical load and in Tables - El for an applied thermal 
load. It is inferred that with increasing gradient index, the non-dimensional 
displacements and stresses increases, whereas, the displacements and the 
stresses decreases with increasing core thickness. This is attributed to the 
change in the fiexural stiffness of the plates due to the increase in metallic and 
ceramic volume fractions. Furthermore, the displacements/stresses predicted 
at the neutral surface using HSDTll and HSDT13 are somewhat different 
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from those of HSDT9 and FSDT5. However, the latter models cannot predict 
the through thickness variation of displacements/stresses. 
Through the thickness variation of displacements and stresses for 1-2-1 Type 
A FGM sandwich plate is shown in Figures [3] - H] for different models pro- 
posed in this study for mechanical and thermal load cases, respectively. For 
mechanical loading case, the higher-order (HSDT13 and HSDTll) and lower- 
order models (HSDT9 and FSDT5) yield almost similar results. The main 
difference among the higher-order elements is in accounting for the slope dis- 
continuity in the in-plane displacements through the thickness, whereas, it 
is the inclusion of variation up to cubic in the in-plane displacements among 
lower-order elements. Since the material properties are smoothly varying 
through the thickness in FGM plates, there is no variation in the evaluated 
results seen among the higher-order and lower-order models. However, the 
performance of these types of elements are significantly different for thermal 
case (see Figure H]), in particular, stress variation and this can be attributed 
to the variation of the coefficient of thermal expansion in the thickness di- 
rection. This effect is captured due to the presence of linear/quadratic terms 
in lateral deflection. 
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Table 2: The effective material properties are based on the Mori-Tanaka 
homogenization scheme. The results are based on the HSDTll element and 
the gradient index, n = 1. 

(a) Comparison of displacements, stresses for the Al/SiC functionally 
graded square plate when subjected to a mechanical load. The nondi- 
mensionalized physical quantities are: u = 100Emu/{qohS^),w — 

100ErnW/{qohS^),{axx,'^xy) = Wia^x, CTxy) / (qoS^) ,axz = lOaxz/{qoS) 

a/h Mesh Varaible 



U W O'xx ^ xy ^ xz 

4x4 -2.9305 2.5407 2.7905 -1.6545 2.2783 

6x6 -2.9153 2.5512 2.7697 -1.5960 2.2985 

5 8x8 -2.9129 2.5535 2.7549 -1.5783 2.3016 

16 X 16 -2.9129 2.5535 2.7549 -1.5783 2.3016 

Ref. [2Tj -2.9129 2.5559 2.7562 -1.5600 2.3100 

4x4 -2.8972 2.1169 2.6564 -1.6470 2.3253 

6x6 -2.8972 2.1148 2.5762 -1.5880 2.3287 

40 8x8 -2.8967 2.1152 2.5494 -1.5704 2.3286 

16 X 16 -2.8967 2.1152 2.5494 -1.5704 2.3286 

Ref. [21] -2.8984 2.1163 2.6093 -1.5522 2.3281 

(b) Comparison of displacements, stresses for the Al/SiC functionally 
graded square plate when subjected to a thermal load. The nondi- 
mensionalized physical quantities are: u — lOOu/{hamToS),w = 

100w/{hamToS'^), {cTxx.^xy, CFxz) = lQ{crxx, CFxy , axz) / {EynO-TnTo) 



a/h Mesh Varaible 



40 



U W^ Oxx CFxy O'xz 

4x4 -1.3176 4.4726 -1.2966 -6.6678 2.0430 

6x6 -1.3232 4.4946 -1.6720 -7.2500 2.3848 

8x8 -1.3240 4.4499 -1.7445 -7.1786 2.5352 

16 x 16 -1.3240 4.4499 -1.7445 -7.1786 2.5352 

4x4 -1.2802 3.3660 -1.6954 -7.2616 2.3157 

6x6 -1.2868 3.3795 -2.0216 -7.0540 2.6674 

8x8 -1.2878 3.3818 -2.1340 -6.9832 2.8544 

16 X 16 -1.2878 3.3818 -2.1340 -6.9832 2.8544 
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Table 3: Deflections and stresses for a simply supported square FGM sandwich plates with homogeneous 
core (Type A), with gradient index, n = 0.5, subjected to a sinusoidally distributed load (mechanical load). 
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O'xx 


f^xy 


O'xz 
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O'xz 




HSDT13 


0.01827 


0.01257 


-0.05962 


0.03131 


0.26308 


0.01677 


0.01158 


-0.05469 


0.02872 


0.26010 
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HSDTll 


0.01827 


0.01257 


-0.05964 


0.03131 


0.26321 


0.01677 


0.01158 


-0.05471 


0.02872 


0.26010 


HSDT9 


0.01867 


0.01342 


-0.05942 


0.03199 


0.26284 


0.01712 


0.01236 


-0.05447 


0.02933 


0.25980 




FSDT5 


0.01813 


0.01358 


-0.05769 


0.03106 


0.26449 


0.01656 


0.01247 


-0.05269 


0.02837 


0.26158 




HSDT13 


0.01818 


0.01181 


-0.05822 


0.03114 


0.26438 


0.01662 


0.01081 


-0.05323 


0.02847 


0.26150 
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HSDTll 


0.01817 


0.01181 


-0.05823 


0.03114 


0.26447 


0.01662 


0.01081 


-0.05324 


0.02847 


0.26150 


HSDT9 


0.01827 


0.01201 


-0.05813 


0.03130 


0.26413 


0.01670 


0.01099 


-0.05313 


0.02861 


0.26119 




FSDT5 


0.01813 


0.01205 


-0.05769 


0.03106 


0.26450 


0.01656 


0.01102 


-0.05269 


0.02837 


0.26159 




HSDT13 


0.01813 


0.01154 


-0.05771 


0.03107 


0.26548 


0.01655 


0.01054 


-0.05270 


0.02838 


0.26255 
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HSDTll 


0.01813 


0.01154 


-0.05771 


0.03107 


0.26548 


0.01655 


0.01054 


-0.05270 


0.02838 


0.26255 


HSDT9 


0.01813 


0.01154 


-0.05769 


0.03107 


0.26515 


0.01655 


0.01054 


-0.05268 


0.02838 


0.26222 




FSDT5 


0.01813 


0.01155 


-0.05768 


0.03107 


0.26514 


0.01655 


0.01054 


-0.05268 


0.02837 


0.26221 



Table 4: Deflections and stresses for a simply supported square FGM sandwich plates with homogeneous 
core (Type A), with gradient index, n = 5, subjected to a sinusoidally distributed load (mechanical load). 



--I 



a/h 


Element 






1-1-1 










1-2-1 








Type 


u 


w 


'^xx 


(^xy 


'^xz 


u 


w 


^xx 


(^xy 


'^xz 




HSDT13 


0.04232 


0.02828 


-0.13876 


0.07250 


0.31653 


0.03233 


0.02151 


-0.10626 


0.05538 


0.31370 


5 


HSDTll 


0.04230 


0.02827 


-0.13871 


0.07247 


0.31626 


0.03232 


0.02151 


-0.10626 


0.05537 


0.31360 


HSDT9 


0.04362 


0.03024 


-0.13883 


0.07474 


0.31516 


0.03332 


0.02303 


-0.10602 


0.05708 


0.31277 




FSDT5 


0.04353 


0.03077 


-0.13852 


0.07457 


0.31534 


0.03277 


0.02340 


-0.10430 


0.05615 


0.31376 




HSDT13 


0.04323 


0.02785 


-0.13860 


0.07407 


0.31601 


0.03267 


0.02102 


-0.10483 


0.05598 


0.31413 


10 


HSDTll 


0.04323 


0.02784 


-0.13860 


0.07406 


0.31599 


0.03267 


0.02102 


-0.10483 


0.05597 


0.31409 


HSDT9 


0.04355 


0.02834 


-0.13859 


0.07462 


0.31534 


0.03291 


0.02140 


-0.10473 


0.05639 


0.31356 




FSDT5 


0.04353 


0.02847 


-0.13852 


0.07458 


0.31537 


0.03277 


0.02150 


-0.10430 


0.05615 


0.31378 




HSDT13 


0.04351 


0.02771 


-0.13853 


0.07460 


0.31667 


0.03277 


0.02086 


-0.10432 


0.05617 


0.31504 


100 


HSDTll 


0.04351 


0.02771 


-0.13853 


0.07460 


0.31667 


0.03277 


0.02086 


-0.10432 


0.05617 


0.31504 


HSDT9 


0.04352 


0.02771 


-0.13849 


0.07460 


0.31626 


0.03277 


0.02087 


-0.10428 


0.05618 


0.31464 




FSDT5 


0.04352 


0.02771 


-0.13849 


0.07460 


0.31623 


0.03277 


0.02087 


-0.10428 


0.05617 


0.31461 
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Figure 3: Displacements and stresses through the thickness for the square 
plates with simply supported edges for 1-2-1 Type A FGM plate with gradient 
index, n = 1, a/h = 5, subjected to a uniformly distributed mechanical 
loading. 
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Table 5: Deflections and stresses for a simply supported square FGM sandwich plates with homogeneous 
core (Type A), with gradient index, n = 0.5, subjected to a sinusoidally distributed load (thermal loading). 



a/h 


Element 
Type 






1-1-1 










1-2-1 








u 


w 


C^xx 


Cxj/ 


^xz 


u 


w 


^xx 


<7xy 


O^xz 




HSDT13 


0.14143 


0.09551 


1.13358 


0.34643 


-0.01628 


0.13275 


0.08966 


1.21775 


0.32518 


-0.00813 


5 


HSDTll 


0.14148 


0.09551 


1.13318 


0.34657 


-0.01637 


0.13280 


0.08965 


1.21740 


0.32529 


-0.00826 


HSDT9 


0.13894 


0.08825 


0.79687 


0.34030 


-0.01026 


0.13038 


0.08282 


0.83580 


0.31933 


-0.00732 




FSDT5 


0.13821 


0.08799 


0.80020 


0.33845 


-0.01046 


0.12962 


0.08252 


0.83924 


0.31742 


-0.00755 




HSDT13 


0.13902 


0.08987 


1.14995 


0.34052 


0.02687 


0.13041 


0.08431 


1.23367 


0.31943 


0.01808 


10 


HSDTll 


0.13903 


0.08987 


1.14982 


0.34055 


0.02692 


0.13042 


0.08431 


1.23355 


0.31946 


0.01811 


HSDT9 


0.13839 


0.08806 


0.79938 


0.33895 


0.01389 


0.12981 


0.08260 


0.83839 


0.31793 


0.00891 




FSDT5 


0.13821 


0.08799 


0.80020 


0.33846 


0.01397 


0.12962 


0.08252 


0.83924 


0.31742 


0.00896 




HSDT13 


0.13821 


0.08801 


1.15545 


0.33851 


0.00275 


0.12962 


0.08254 


1.23900 


0.31747 


0.00185 


100 


HSDTll 


0.13821 


0.08801 


1.15545 


0.33851 


0.00275 


0.12962 


0.08254 


1.23900 


0.31747 


0.00185 


HSDT9 


0.13821 


0.08799 


0.80022 


0.33850 


0.00142 


0.12962 


0.08253 


0.83927 


0.31747 


0.00092 




FSDT5 


0.13820 


0.08799 


0.80023 


0.33849 


0.00142 


0.12962 


0.08252 


0.83927 


0.31746 


0.00092 



to 

o 



Table 6: Deflections and stresses for a simply supported square FGM sandwich plates with homogeneous 
core (Type A), with gradient index, n = 5, subjected to a sinusoidally distributed load (thermal loading). 



a/h 


Element 
Type 






1-1-1 










1-2-1 








u 


w 


C^xx 


Cxj/ 


^xz 


u 


w 


^xx 


<7xy 


C^xz 




HSDT13 


0.18232 


0.12238 


0.73973 


0.44652 


-0.02475 


0.15938 


0.10641 


0.96614 


0.39045 


-0.05416 


5 


HSDTll 


0.18240 


0.12247 


0.73915 


0.44679 


-0.02497 


0.15942 


0.10644 


0.96585 


0.39059 


-0.05439 


HSDT9 


0.17878 


0.11317 


0.61575 


0.43787 


-0.01539 


0.15555 


0.09815 


0.72136 


0.38102 


-0.03237 




FSDT5 


0.17771 


0.11314 


0.62060 


0.43519 


-0.01570 


0.15382 


0.09793 


0.72923 


0.37668 


-0.03299 




HSDT13 


0.17888 


0.11544 


0.76316 


0.43812 


0.02160 


0.15522 


0.10004 


0.99415 


0.38025 


0.01437 


10 


HSDTll 


0.17889 


0.11547 


0.76301 


0.43819 


0.02166 


0.15523 


0.10005 


0.99408 


0.38029 


0.01442 


HSDT9 


0.17798 


0.11315 


0.61940 


0.43590 


0.01152 


0.15425 


0.09798 


0.72728 


0.37783 


0.00732 




FSDT5 


0.17771 


0.11314 


0.62060 


0.43519 


0.01159 


0.15382 


0.09793 


0.72924 


0.37668 


0.00736 




HSDT13 


0.17771 


0.11317 


0.77116 


0.43527 


0.00221 


0.15382 


0.09795 


1.00359 


0.37676 


0.00147 


100 


HSDTll 


0.17771 


0.11317 


0.77116 


0.43527 


0.00221 


0.15382 


0.09795 


1.00359 


0.37676 


0.00147 


HSDT9 


0.17771 


0.11314 


0.62064 


0.43525 


0.00117 


0.15381 


0.09793 


0.72926 


0.37674 


0.00075 




FSDT5 


0.17770 


0.11314 


0.62064 


0.43524 


0.00117 


0.15381 


0.09793 


0.72927 


0.37672 


0.00075 
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Figure 4: Displacements and stresses through the thickness for the square 
plates with simply supported edges for 1-2-1 Type A FGM plate with gradient 
index, n = 1, a/h = 5, subjected to a thermal loading. 



4.2 Free flexural vibrations 

Next, the free vibration characteristics of FGM sandwich plate is numeri- 
cally studied. In all cases, we present the non-dimensionalized free flexural 
frequency defined as: 



Q 



uja 




(22) 



where uj is the natural frequency, Po = 1 kg/m^ and Eo=l GPa. Before pro- 
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ceeding with the detailed study on the effect of gradient index and the type 
of sandwich FGM plates on the natural frequencies, the formulation devel- 
oped is validated against available three-dimensional elasticity solutions [11]. 
Based on a progressive refinement, a 8x8 mesh is found to be sufficient to 
model the sandwich FGM plate. Table [7] gives a comparison of the first six 
computed frequencies for a simply supported square Type A FGM sandwich 
plate. It can be seen that the numerical results from the present formulation 
are found to be in very good agreement with the existing solutions. 
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Table 7: Convergence of flexural vibration frequency parameters Q of square 
2-1-2 FGM plates of Type A. 

n a/h Mesh Frequency 



10 



10 



ill * '2 ''3 ''4 * '5 ' '' 



6 



4x4 1.2297 2.6874 2.6876 2.8149 2.8271 4.1312 

6 x6 1.2294 2.6869 2.6869 2.8032 3.2551 4.1553 

8 x8 1.2293 2.6868 2.6868 2.8009 2.8345 4.1568 

16 xl6 1.2293 2.6868 2.6868 2.8009 2.8345 4.1568 

4x4 1.3025 3.1822 3.1822 4.9167 6.2782 7.7331 

6 x6 1.3020 3.1640 3.1640 4.9209 6.0938 7.6898 

8 x8 1.3019 3.1606 3.1606 4.9188 6.0586 7.7675 

16 xl6 1.3019 3.1606 3.1606 4.9188 6.0586 7.7675 

Li et al, [llj 1.3018 3.1588 3.1588 4.9166 6.0405 



4x4 0.8958 2.0729 2.0729 2.2067 2.2225 3.0703 

6x6 0.8955 2.0636 2.0637 2.2063 2.5395 3.0840 
5 

8x8 0.8955 2.0618 2.0619 2.2062 2.4973 3.0845 

16x16 0.8955 2.0618 2.0619 2.2062 2.4973 3.0845 

10 4x4 0.9423 2.3109 2.3109 3.5841 4.5875 5.6835 

6x6 0.9419 2.2974 2.2974 3.5851 4.4488 5.6311 

8x8 0.9418 2.2948 2.2948 3.5832 4.4225 5.6183 

16x16 0.9418 2.2948 2.2948 3.5832 4.4225 5.6183 

Li et al, fTT] 0.9404 2.2862 2.2862 3.5647 4.3844 
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Table 8: Fundamental frequency parameters Q of simply supported square FGM sandwich plates with 
homogeneous core (Type A). 

a/h Element 1-1-1 1-2-1 2-2-1 

Type 0.5 1 5 0.5 1 5 0.5 1 



5 



to 



HSDT13 1.6774 1.4219 1.2778 0.9986 

HSDTll 1.6774 1.4219 1.2778 0.9988 

HSDT9 1.6774 1.4152 1.2714 0.9937 

FSDT 1.6689 1.4076 1.2628 0.9860 

Li et al, [H] 1.6771 1.4218 1.2777 0.9980 



1.4696 1.3536 1.1192 

1.4696 1.3537 1.1193 

1.4626 1.3468 1.1131 

1.4565 1.3398 1.1053 

1.4694 1.3534 1.1190 



1.4455 1.3144 1.0565 

1.4455 1.3144 1.0566 

1.4387 1.3078 1.0510 

1.4320 1.3002 1.0444 

1.4454 1.3143 1.0561 



HSDT13 
HSDTll 
10 HSDT9 
FSDT 



1.8269 1.5214 1.3553 1.0455 

1.8269 1.5214 1.3553 1.0456 

1.8245 1.5193 1.3553 1.0441 

1.8242 1.5168 1.3506 1.0418 



Li et al, [11] 1-8268 1.5213 1.3552 1.0453 



1.5768 1.4415 1.1757 

1.5768 1.4415 1.1758 

1.5746 1.4394 1.1740 

1.5726 1.4371 1.1715 

1.5767 1.4414 1.1757 



1.5494 1.3977 1.1100 

1.5494 1.3977 1.1100 

1.5472 1.3957 1.1084 

1.5451 1.3932 1.1064 

1.5493 1.3976 1.1098 



HSDT13 
HSDTll 
100 HSDT9 
FSDT 



1.8884 1.5605 1.3852 1.0631 

1.8884 1.5605 1.3852 1.0631 

1.8883 1.5605 1.3851 1.0631 

1.8883 1.5605 1.3851 1.0631 



Li et al, [H] 1.8883 1.5605 1.3851 1.0631 



1.6192 1.4756 1.1970 

1.6192 1.4756 1.1970 

1.6192 1.4756 1.1970 

1.6192 1.4756 1.1970 

1.6192 1.4756 1.1970 



1.5904 1.4300 1.1303 

1.5904 1.4300 1.1303 

1.5904 1.4300 1.1302 

1.5904 1.4299 1.1302 

1.5903 1.4299 1.1302 



to 



Table 9: Fundamental frequency parameters Q of simply supported square FGM sandwich plates with FGM 
core (Type B). 

a/h Element 1-1-1 1-2-1 2-2-1 

Type 0.5 1 5 0.5 15 0.5 1 5 



10 



100 



HSDT13 1.0893 1.1511 1.1701 1.2162 

HSDTll 1.1078 1.1512 1.1705 1.2184 

HSDT9 1.1021 1.1449 1.1639 1.2113 

FSDT 1.1263 1.1503 1.1642 1.2050 



1.1663 1.1952 1.2712 

1.1664 1.1953 1.2718 
1.1597 1.1884 1.2644 
1.1660 1.1880 1.2567 



1.2031 1.2421 1.3312 

1.2034 1.2422 1.3326 

1.1965 1.2350 1.3249 

1.1950 1.2299 1.3173 



HSDT13 1.2087 1.2392 1.2524 1.2935 

HSDTll 1.2156 1.2392 1.2525 1.2942 

HSDT9 1.2138 1.2373 1.2506 1.2921 

FSDT 1.2225 1.2394 1.2509 1.2903 



1.2598 1.2806 1.3513 

1.2598 1.2806 1.3515 

1.2578 1.2785 1.3492 

1.2601 1.2786 1.3469 



1.2865 1.3238 1.4180 

1.2866 1.3238 1.4184 
1.2846 1.3216 1.4161 
1.2842 1.3201 1.4136 



HSDT13 1.2616 1.2751 1.2854 1.3239 

HSDTll 1.2617 1.2751 1.2854 1.3239 

HSDT9 1.2617 1.2751 1.2854 1.3239 

FSDT 1.2618 1.2751 1.2854 1.3239 



1.2981 1.3148 1.3825 

1.2981 1.3148 1.3825 

1.2981 1.3148 1.3825 

1.2981 1.3148 1.3825 



1.3198 1.3559 1.4519 

1.3198 1.3559 1.4519 

1.3198 1.3559 1.4519 

1.3198 1.3559 1.4518 



Tables [H] and M gives the fundamental frequencies for Type A and Type B 
sandwich FGM plate for three different thickness ratios and for different core 
thickness, respectively. It can be seen that the natural frequencies decreases 
with increasing gradient index for Type A and the natural frequencies de- 
creases with decreasing gradient index for Type B FGM sandwich plate. 
This can be attributed to the decrease in the material rigidity. In the case 
of Type A, as the gradient index increases, the metallic volume fraction in- 
creases, thus decreasing the overall stiffness of the plate. In the case of Type 
B, decreasing the gradient index, decreases the metallic volume fraction and 
increases the material rigidity. 

Figures - El shows the relative displacements through the thickness of the 
six modes using HSDT13. The plate is simply supported square 1-2-1 FGM 
plate with gradient index n = 2 and a/h = 5. The displacements {u, v, w) are 
plotted along the lines (a/2, 6/2, z) and (a/4, 6/4, z), where —h/2 < z < h/2. 
It can be seen that, in flexural mode 1, the transverse displacement w is 
not uniform at the chosen locations, thus, exhibiting the existence of normal 
stresses in the thickness direction. In flexural modes 2-6, the deflected 
shape retains the thickness at (a/2, 6/2, z), whereas the thickness of the plate 
is compressed at the other location, i.e., at (a/4,6/4,2;). It can be seen that 
the variation of the in-plane displacements {u,v), in general, linear or non- 
linear in some higher-order models, irrespective of the types of modes. 

5 Conclusion 

FGM sandwich plate bending and free vibration analyses are carried out con- 
sidering various parameters such as the material gradient index, the sandwich 
type and the thickness ratio. Different plate models are employed in predict- 
ing the physical behaviour and their through thickness variations in the plate. 
The accuracy and the effectiveness of the higher-order models (HSDT13 and 
HSDTll) over the lower-order theories have been demonstrated considering 
problems for which analytical/numerical results are available in the litera- 
ture. It may be concluded that the IISDT13 model predict accurate results 
for any type of sandwich construction whereas other type of models depend 
on the type of sandwich plate and loading situations. 



26 









1 


u 


, 






w 














1 






' 












1 
















. 


I 



1r 
0.8 
0.6 
0.4 

:■ "-2 

(0 

* 

c - 

o 

j= -0.2 

-0.4 



0.4 0.6 

Displacements 

(a) mode 1 









/ 




u 






— e— V 




/ 




w 




/ 










* 




/ 




/ 










/ 


- 


















/ 




f 






- 


/ 













Displacements 



(c) mode 3 







i J y/" 




u 

w 


/ X 




/ 


A 



Displacements 

(e) mode 5 



^% 




• 


u 


-e- 


— e— V 


\ 






\ 




- 






\ 


s 








\ 












' N 














i ■». 




V 












\ 











-0.15 -0.1 -0.05 0.05 0.1 

Displacements 



(b) mode 2 



1 






\ ' y^ 


0.8 








0.6 




w 


\ I y^ 


0.4 






\\ / 


0.2 


\/^ 





ji 


0.2 


y\ 


0.4 


/^ ;\ 


0.6 


y^ l\ 


0.8 

-1 


y^ \ 



Displacements 

(d) mode 4 











^^^ 


u 


0.8 


^^ ; 


w 




^^h ' 




0.4 


\. ; 


0.2 


\ ; 





\ 


0.2 


;\ 


0.4 


; >. 


0.6 


' ^^ 


0.8 

-1 


Xv 



Displacements 

(f) mode 6 



Figure 5: Deflected shapes Ui{a/2,a/2,z), i = x,y,z oi the six modes for 
the square plates with simply supported edges for 1-2-1 Type A FGM plate 
with gradient index, n = 1 and thickness, a/h = 5. 
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Figure 6: Deflected shapes 'Ui(a/4,a/4, 2;), i = x,y,z oi the six modes for 
the square plates with simply supported edges for 1-2-1 Type A FGM plate 
with gradient index, n = 1 and thickness, a/h = 5. 
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